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Abstract 
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I. INTRODUCTION 



Although the standard model (SM) of particle physics consists of three fermion genera- 
tions, the number of these generations is not fixed by theory. Asymptotic freedom in quan- 
tum chromo dynamics (QCD) limits the number of generations to fewer than nine. Neutrino 
counting, based on data at the intermediate vector boson Z, shows that the number of gen- 
erations having light neutrinos (mj, ^ vnzjl ) is equal to 3. However, neutrino oscillations 
suggest a new mass scale that is beyond that of the SM, and the possibility of additional 
heavier neutrinos is open. In addition, it has recently been emphasized that the electroweak 
oblique parameters do not exclude a fourth generation of chiral fermions [lijsl. In the era of 
the Large Hadron Collider (LHC), the search for heavy quarks and leptons beyond those of 
the SM should be kept in mind. 

Observation of new heavy fermions requires knowledge of their expected production cross 
sections and decay properties. In this work, we present inclusive cross sections of fourth 
generation fermion production at three LHC energies, either in pairs or singly, calculated to 
next-to-leading (NLO) order accuracy in QCD. Leading order Feynman diagrams for new 
heavy fermion production are shown for illustrative purposes in Fig. [H We do not show the 
full set of NLO diagrams used in our calculations. Theoretical uncertainties associated with 
higher-order perturbative contributions and with the choice of parton distribution functions 
(PDFs) are specified. We show that the cross sections can be fitted with a rather simple 
analytic formula, and we present the functional form and values of the parameters for the 




FIG. 1: Representative leading order Feynman diagrams for new heavy fermion production. Up- 
percase (lowercase) letters denote new heavy (SM) fermions. The usual gauge bosons are denoted 
byg,PF,Z,7. 
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heavy fermion production cross sections, and their dependence on the mass of the heavy 
fermion. 



II. HEAVY QUARK PAIR PRODUCTION IN QCD 



Heavy quarks that carry the color charge exist in many mode 
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ude models with a sequential fourth generation , or with a vectorli 



s of new physics (NP). 
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11| . or with a b prime (B) [12], or with a heavy top quark partner 
According to their physics motivation, the NP models can be further categorized into two 
classes. In the first group, heavy quarks are introduced to explain discrepancies between 
data and SM predictions. In the second, heavy quarks are essential to solve theoretical 
problems. For example, in Little Higgs models jisl, the heavy top partner cancels the 
quadratic divergence in the quantum corrections to the Higgs boson mass associated with 



the SM top quark loop, 
been studied extensively 



he phenomenology of heavy quark identification at the LHC has 
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25|. 



Similar to the top quark in the SM, heavy quark Q can be pair produced via the strong 
interaction processes, 

qq ^ g ^ QQ and gg QQ, 

where g denotes a gluon. In this work we assume the coupling g-Q-Q is the same as the 
gluon-quark interaction in the SM. The cross section for a different coupling strength can 
be obtained easily from our results by rescaling. In Fig. [2t^a) we plot our calculated NLO 
inclusive cross sections for QQ production at the LHC for three choices of the proton-proton 
center mass (cm.) energy. ^ The square symbol denotes the results of our exact NLO 
calculation whereas the curve is drawn from the simple fitting functions discussed below. As 
shown by the ratios presented in Fig. EJ^b), an increase in the cm. energy from 7 to 14TeV 
can enhance the total cross section markedly, depending on the value of rriQ. The values of 
the standard model parameters used in our calculations are found in Appendix A. 

There are uncertainties in the predicted NLO cross sections associated with the choices 
of the renormalization scale (fiR) and the factorization scale {^f) and with the choice of 
the PDFs. We elect to use the CTEQ6.6M PDFs H for our evaluations. Other sets of 



Our numerical code for the NLO QQ inclusive cross section is based on the analytic results in Refs. 0-0 
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FIG. 2: (a) NLO QQ pair production cross sections as a function of the quark mass mq at the 
LHC. The square symbols denote the results of the exact calculation, while the curves present 
results of our phenomenological fit. (b) Ratios of the cross sections at different cm. energies, 
(c-e) Theoretical uncertainties of the QQ production cross section. The band denotes the PDF 
uncertainties, while the black solid (red dashed) curve denotes the scale dependence obtained by 
varying the renormalization scale by a factor 2 about the central value = t^q- 



PDFs are available, notably the MSTW 3l| and NNPDF 32| sets, and one might use these 
sets in addition to the CTEQ6.6M set to compute a range of predictions. All of these PDF 
sets include an analysis of the uncertainties in the PDF determinations associated both with 
uncertainties in the data used in the global fits and with the choices of theoretical expressions 



used to fit the data. A comparison may be found in Ref. [3l||. In this paper, we compute 
the PDF uncertainty (at 90% CL.) from the master formula in Eq. (2.5) in Ref. 30|], using 
all 44 sets of the CTEQ6.6M package. For heavy quark production at the large masses we 
are considering, the PDF uncertainties arise primarily from the valence quark PDF where 
different PDF sets tend to agree better than they do for gluon distributions. 

In Fig. [2I^c)-12]^e) we display the uncertainties associated with the choices of PDFs and the 
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renormalization/factorization scale. ^ The uncertainties are portrayed relative to the cross 
section with the best-fit PDF by the bands in Fig. [2l The PDF uncertainties are relatively 
large because the mass of the heavy quark requires that the PDFs be sampled in regions 
where they are relatively unconstrained. The PDF uncertainties decrease with increasing 
cm. energy. 

The uncertainties in the NLO cross section associated with the renormalization scale 
{fiji) and factorization scale (fip) are shown in Fig. |2j These uncertainties can be considered 
as an estimate of the size of unknown higher-order contributions. In this study, we set 
A* = f^R = f^F and vary it around the central value of fiQ = mq, where rriQ is the mass of the 
heavy quark. Typically, a factor of 2 is used as a rule of thumb, and we display curves with 
/i = 2/io and fi = ^o/2. In Figs. I21c)-12]^e) we plot 1 -|- o"(/ij)/cr(/io) as a function of mg. The 
cross sections vary between about —10% for fi = 2yUo and +10% ~ 15% for fi = fio/2. We 
note that the scale dependence at the three energies is insensitive to rriQ. It is comparable 
to the PDF uncertainties in the region of relatively small mq but is much smaller than the 
PDF uncertainties in the region of large rriQ. 

We find that our calculated cross sections can be fitted well by a simple three parameter 
analytic expression: 

= ^W(^)-' + BW + c(,)(^). (1) 

where the units of a and rriQ are picobarn (pb) and TeV. The variation with respect to the 
reference cross section (corresponding to fi = mq) gives the NLO scale dependence A^(mQ). 
The values of A^, and coefficients A, B, C are listed in Table [H Our results are presented 
in the form 

- = -(/^o)i:S^;/?, (2) 

where o"(/io) is our prediction at the scale /iq. The cross sections at the scale /io/2 and 2/io 
read as cr(/io/2) = o-(/Uo) + A(j(/io/2) and cr(2/io) = cr(|Uo) + A(j(2yUo), respectively. 

As is evident in Fig. [2]^a), the fitting function and parameters in Table |T] provide an 
excellent representation of the calculated cross sections over the range of heavy quark masses 

^ For the scale dependence, Acr/cro is defined as — (t(/io))/o'(mo), while for the PDF dependence 

Act is defined as Act = 1/4 ^^"('''i^ ~ '^lY^ where where denotes the upper (lower) cross section for 
the ith eigenvector of the PDF sets and n — 22 for CTEQ6.6M PDF sets. The PDF uncertainties are 
then plotted as 1 ± /^a/aQ where ctq stands for the cross section of the averaged PDF set. This definition 
applies to all other figures with cross section uncertainties. 
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FIG. 3: Fractions of the qq and gg initial state contributions to the NLO total cross sections for 
QQ pair production at the LHC. 

shown in the figure. We caution that the expression should not be used outside this range. 
One could extend the region of applicability of a fit to lower mass, say down to the top 
quark mass, by simply increasing the number of terms in the polynomial expression. In 
Appendix |B] we extend our fit of the cross sections for heavy quark pair production to cover 
the wider mass range, 250 to 700 GeV. We also present a tabulation showing our exact NLO 
calculation and our fitted results. The numerical comparison shows that our fitted values 
agree with the exact calculations well within 1%. Interested readers may contact the authors 
for accurate fitting functions that cover the mass range from the top quark mass to 4 TeV. 

One reason more terms may be needed in a polynomial expression is that processes with 
initial state gluons play an important role in pair production of relatively light quarks. As 
shown in Fig. |3l at ^/s = 7 TeV, the gluon initial state contribution dominates over the 
quark initial state when rriQ < 550 GeV; at ^/s = 10 TeV when rriQ < 750 GeV; while at 
^/s = 14 TeV when rriQ < 1 TeV. For relatively light quark production, say rriQ < 300 GeV, 
and terms {y = mg/TeV) are required to fit the total cross section. In Appendix [Ct 
we explore the apparent process independence of the remarkably simple analytic expression 
Eq. [1] and discuss its range of applicability in the heavy fermion mass. 

In view of the interest in the top quark as the heaviest of the known quarks, we provide 
NLO predictions of the top quark pair production cross sections in Table Hi at three LHC 
energies for a few values of m^. Our result at 14 TeV is consistent with Ref. 33|. 
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TABLE I: Fitting parameters {A, B, C) of the parametric formula [Eq. [T] for the NLO total 
cross sections (pb) for heavy quark production at the LHC. The central values are given for the 
scale choice = fip = fJ-o = rnq. The superscript and subscript denote the deviations in the 
coefficients from the central value when the scale is /io/2 (superscript) or 2^o (subscript). We use 
the CTEQ6.6M PDFs. 



Process 


Parameter 


= 7TeV 


^ = 10 TeV 


^= 14 TeV 


qq/99 - 


-+ QQ 


A 
B 
C 


0.84569l°:S 

i.48bt)t)_Q 02244 
-o-o'U48_o.08488 


i./DDDi_o.01912 

-0 391 94+0-00546 
u.oyiy4_o 06751 

c; 091/19+O.O4314 
O.OZi4Z_o. 03726 


9 nQQQQ+0.02345 
Z.UOOOO_o 01292 

-0 45930+°-°3295 

U.40you_o 08482 

4 4^^8^^^+0-01784 
-4.4£)»t)d__o.oi642 


qb q'T 


A 
B 
C 


n 79090+0.00130 

U. (zozo_o 01262 

1 fifi"?! 8~0 03563 
l.OUOiO_,_o 00630 

q qoqo/i +0.05251 
0.00004_o 04100 


72901-0-00419 

U./zjUi_,_0.00981 

9 9Sqf;4-0-01591 
Z.ZOOO4_,_0 00541 

9 c;7c;nQ+0.03980 
Z.orouo_o 03535 


0.69166;°:™ 

9 04970-0.01395 
z.y^z/y_o 01098 

-2.06090j:°:°i}^ 


qb q'f 


A 
B 
C 


n 777/19+0.00580 
U. i 1 '4Z_Q 00158 

n Sfi99Q-0-04447 
u.ouzzy_^0 03096 

q r7qc-q+0. 05666 
0.d(0<JO_o 05169 


n 7«fi1 7+0.00133 
U-'^OJ-' +0.00520 

1 c;qqQ7-0.03343 
i.OOOOr_^0.01617 

9 797Q1 +0.04759 
z. (Z(yi_o 03900 


n 7c^S9Q-0.00266 
U. 1 Oozy_^o 00743 

9 99Qoq-0. 01720 
Z.zzyoo_,_0 01447 

9 1 8987+0.03754 
Z.10ZO(_o 03617 


qq' - 


> Tb 


A 
B 
C 


1 071 nq-0. 00541 
1.0/ 1JO_^0. 01703 

J.44D/i„o.05685 

O.ZZyOD_o.oi643 


1 41641+0-00152 
i.41D4i_,_0 00749 

q 4c;Qnf;+0-01625 
0.400UO_o 02949 

Q 01 qni +0.02269 
O.OiOU±_o,oi891 


1 4692^-°-°°'^^3 

1.4Dyzo_^0 00596 

T qqc:qQ+0.02351 
O.OOOOy_o.02316 

-3.00593t°:°}i8 


qq' - 


> fb 


A 
B 
C 


1-79403;0:S 

-5.40590t°:S 
-4.81694l°:S 


1-77007;°:™ 
-5.03094tO:S 

q f;7Qq-| +0.01489 
O.U(001_o 01460 


1 fi8Q4n-0-00650 
1.68940^0.00815 

-4 44604+0-02719 

4.44DU4_o 02747 

nS8n'?+0-00993 
-3.U8803_o.oio88 


qq' - 


' TB 


A 
B 
C 


u.uyouo_,_o. 00202 

1 nSQI ^^+0-03084 
-i.U»yit)_o.03419 

1 n 4709+0.04744 

-iU.4/82_o.o4719 


0.98804l°:S 

9 9i9trc-0.01419 
Z.Z1ZOD_0.00969 

A ««oin+0.04337 
-4.66219_o.o3293 


1 0701 0-0 00263 
i.UmiU_,_0. 00304 

-2 20446+°-°^^^^ 

Z.ZU44D_o.oi452 

7c;9i 7+0.01638 
0. rOZi r_o 01956 


qq' - 


' TB 


A 
B 
C 


n fif^7'^8+0-01149 
U.btiM»_0 00549 

-4 021 7'S"0-0^027 
4.UZl/D_o.o0550 

-9.05880lo:o6865 


1 1^^00-0-00182 
i.iOOUU_,_0.00280 

r, nnoqi +0.02561 
0.uyoOi_o.02683 

-6.08544t°:S 


1 1 8fifi'^-0 00354 
l.lOUOO_,_0 00487 

A 7qq/|o+0.02658 
-4. M348_o 02800 

-4 91 903+0-01753 
^.yiyuo_o.oi835 
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TABLE II: NLO top quark pair production cross sections (pb) at the LHC with three different 
cm. energies based on the CTEQ6.6M PDFs. The central values are given for the scale choice 
IJ'R = = /^o = "T-t- The superscript and subscript denote the deviations of the coefficients from 
the central value when the scale is //o/2 (superscript) or 2/io (subscript). 



mt(GeV) 



171 
173 
175 
177 



^ = 7TeV y/s = lOTeV y/s = 14TeV 



m.AtTz 396.5l^s^i 877.2l}°?i 



148.21},^^ 374.5^^1^ 



830.9" 



-93.9 
-94.1 



139.6 



-16.5 
-18.1 



354.0 



4-41.6 
-43.4 



787.5+!^j-? 



131.61}^:^ 335.01^?:? 747.2 



91.5 
-1-86.4 



III. HEAVY QUARK PRODUCTION THROUGH THE ELECTROWEAK IN- 
TERACTION 



Heavy quark production occurs also through electrowea 



i interactions. One example is 



the heavy top quark partner (T) in the Little Higgs Model 1^, [isl, which is responsible for 
canceling large quantum corrections to the Higgs boson mass from the SM top quark loop. 
In addition to its pair production through QCD interactions, the T quark can be produced 
singly in association with a light quark through the processes 

qh —7- g'T, and qq' — )■ Th. 

It may also be produced in association with a heavy B quark via the s-channel subprocess 

qq' ^W* ^ TB. 

We consider first single T quark production in association with a SM h quark via an 
s-channel subprocess, illustrated at leading order in Fig. Hl^b). This subprocess involves a 
weak interaction vertex W-T-h and a quark mixing matrix Vrb- For simplicity, we assume 
that the coupling of W-T-h is same as the SM W-t-b coupling and that | Vrb] = 1. Our results 
may be rescaled for other choices. We calculate the Tb and Tb results separately since the 
pp initial state at the LHC is not a CP eigenstate. The renormalization and factorization 
scales are chosen to he fiR = fip = \/ Q\ where Qy denotes the four- momentum of the W 
boson in the propagator. 

The square symbols in Fig. H] show the predicted production cross sections for Tb and Tb 
at the LHC for three cm. energies. The curves are our fit with the parametrization in Eq. [TJ 
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FIG. 4: NLO production cross sections for (a) Tb and (b) Tb at the LHC as a function of m-r, 
and (c-e) PDF and scale uncertainties of the Tb production cross section. The Tb production cross 
section exhibits similar uncertainties. 

We note that this simple parametrization fits single T production via the s-channel process 
very well for the three LHC energies. The fitting parameters are presented in Table HI The 
cross section for single T production is larger than the cross section for single T production 
by almost a factor 2 ~ 3. This is due to the difference in parton densities of the colliding 
protons. While in both cases the antiquark is from the quark sea of one of the incoming 
protons, the probability that it collides with an up quark from the other proton is higher 
than the probability for a collision with a down quark. Since Tb and Tb production exhibit 
similar uncertainties, we present only the results for Tb production in Figs. Hl^c-e). We note 
the scale dependence at the three energies is not very sensitive to m^, and it is much smaller 
than the PDF uncertainties in the large region. An increase in the cm. energy reduces 
the PDF uncertainties in the large rriT region sizably. 

In the "sequential" fourth generation model, the heavy T and B quarks form an isodou- 
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FIG. 5: NLO production cross sections for (a) TB and (b) TB at the LHC as a function of 
niQ = rriT = ms, and (c-e) PDF and scale uncertainties of the TB production cross section. The 
TB production cross section exhibits similar uncertainties. 

blet, and their interaction with the W boson is identical to the W-t-b interaction in the SM 
with the substitution of Vtb ^ Vtb = ^- For simphcity, we take the T and B quarks to be 
degenerate, i.e. mx = ms = mq. The T and B quarks can be produced together via the 
Drell-Yan process (with an s-channel W*). In Fig. [5] we show our calculated cross sections 
and the fitted curves based on the parametrization in Eq. [1] The parametrization works 
well at ^/s = 10 TeV and 14 TeV but less adequately at 7 TeV. The inclusion of more terms 
in the polynomial would improve the fit. The fitting parameters are listed in Table [H 

Last, we consider single T quark production via a t-channel exchange diagram, illustrated 
at leading order in Fig. (U^c). The production cross section is larger than in the s-channel 
case because it is not subject to 1/s suppression. In this calculation we adopt the so-called 



double deep inelastic-scattering scale proposed in Ref. [sj]- Color conservation enforces a 



natural factorization of the scales. The fermion line in Fig. [T]^c) that does not include a heavy 



quark probes a proton with the deep inelastic-scattering scale a/ Qy , which is identical to 
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FIG. 6: Single heavy (a) T production and (c) T production via the t-channel subprocess as a 
function of rriT- The square symbols denote the results of the exact calculation while the curves 
present results of our phenomenological fit. (b) and (d) Ratios of the cross sections at different 
cm. energies, (e-g) Theoretical uncertainties of the single T production cross section. The band 
denotes the PDF uncertainties while the black solid (red dashed) curve denotes the scale dependence 
obtained by varying the renormalization scale by a factor 2 about the central value /io = mT- 

the virtuality of the W boson through NLO. The fermion line that connects to a heavy 
quark sees the deep inelastic-scattering scale for massive quarks of y^Q^H-m^. 

Our simple parametrization also works well for the t-channel results as may be seen in 
Figs. El^a) and El(c). The fitting parameters are found in Table [B Again, an increase in 
the cm. energy from 7 to 14 TeV can enhance the total cross section markedly, depending 
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on the value of m^; see Figs. |6]^b) and[6]^d). The PDF and scale uncertainties are shown 
in Figs. |6t^e)-l6t^g). Since single T and single T production exhibit similar uncertainties, we 
present only the results for single T production. We note the scale dependence at the three 
energies is not sensitive to m^. The cross sections vary between ±2% for rriT ~ 500 GeV to 
±7% for rriT ~ 2 TeV. The scale dependence is much smaller than the PDF uncertainties in 
the large rriT region. 



IV. EXOTIC LEPTON PRODUCTION 



Several variants of the seesaw mechanism have been proposed to explain light neutrino 
masses. Positive signals could be observed at the LHC in the case that seesaw messengers 
exist at the TeV scale or below. Three types of tree-level seesaw models generate light 
neutrino Majorana masses via introduction of (1) a right-handed neutrino singlet, (2) a 
complex scalar triplet A with hypercharge Y = 1, and (3) a lepton triplet E with Y = 0. In 
these seesaw models, neutrinos achieve their masses via a lepton number violating operator 
at the scale A. The scale is not necessarily very high, and it might be around a TeV. Searches 
for these exotic leptons and, if observed, measurement of their properties would verify and 
even distinguish different seesaw mechanisms ^|. For this purpose, an accurate calculation 
including higher-order QCD corrections is in order. 

TABLE III: Coupling strengths of the gauge vertices in heavy lepton pair production for differ- 
ent models, where g denotes the usual weak coupling strength. Note that all the couplings are 
vector like. 

WEN ZEE -lEE ZNN 

Majorana triplet g gc^r e 

Dirac triplet g gc^ e 

Lepton isodoublet -^( h Svf) e ^ 

^/2 Cw 2 2cw 

Charge singlet ... ... e 
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We examine the following three possibilities for exotic lepton production: 

qq^-i/Z E+E-, 
qq^ Z ^ NN, 
qq' ^ EN, 

where E{N) denotes the new heavy charged lepton (heavy neutrino). The interactions of 
gauge bosons and heavy leptons are summarized in Table UTTl for four interesting NP models: 
(1) a Majorana lepton triplet with Y = 0, (2) a Dirac lepton triplet with Y = 0, (3) a lepton 



isodoublet, (4^ 
found in Ref. 



a charge singlet. The motivation for these models and more details can be 



361]. It is worth mentioning the following properties: 



All the interactions given in the Table IIIII are vectorlike. 

There are no ZNN interactions in the Majorana/Dirac lepton triplet model because 
the triplet has zero hypercharge, and NN pairs are not produced. 

A Dirac lepton triplet is forme d by two degenerate Majorana triplets with opposite 
CP parities. As shown in Ref. 



37l |. the heavy lepton fields can be redefined in such 
a way that the Lagrangian is written in terms of two charged leptons E^ , E2 and a 
Dirac neutrino A^. As a consequence of the presence of two charged fermions instead of 
only one, the total heavy lepton production cross section is twice that for a Majorana 
triplet. In this work, however, we consider only one flavor of charged lepton in both 
the Majorana triplet and the Dirac triplet. 

To make our results more useful, rather than focusing on a specific NP model, we consider 
the following effective interaction in our numerical calculation, 

C = -— (EYE + N-f^'N) Z^--^ (E-i^'NW- + N-f^'EW^) + eE-^^'EA^, (3) 

where the symbol g stands for the usual weak coupling strength. Numerical results for the 
models listed in Table IIIII can be easily derived from ours by simply rescaling the coupling 
strengths, except for the process gg — )■ 7/Z — )■ E^E~ where there are interference effects of 
the photon and Z boson contributions. However, since the threshold for two heavy charged 
leptons is much higher than the Z boson mass, the interference effects are subleading. We 
present separately the NLO cross sections for qq ^ j ^ E^E^ and qq ^ Z ^ E^E^; the 
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FIG. 7: NLO production cross sections for (a) E^N and (b) E^N production at the LHC as a 
function of rriL = tue = mi^, and (c-e) PDF and scale uncertainties of E~^N production. The 
E^N production cross section exhibits similar uncertainties. 

couplings involved in the two processes are overall factors. For all the s-channel processes 
listed above, the renormalization and factorization scales are chosen to be = \/ Q\ where 
Qy denotes the four-momentum of the gauge boson in the propagator. For simplicity, we 
take the E and leptons to be degenerate, i.e. rriE = Tn^q = rriL. 

In Figs, d^a) and[71^b) we display the NLO cross sections for the qq' — )■ — ?■ E^N and 
qq' — 7- W~ — E~N processes, respectively. Both processes exhibit similar scale and PDF 
uncertainties as shown in Figs. [7] (c)-me). The fitting parameters are given in Table HVl 

We also plot the NLO cross sections for the gg — t- 7 — ?■ E^E~ and the qq ^ Z ^ E^E~ 
processes in Fig. [HI Both processes exhibit similar scale and PDF uncertainties as shown in 
Fig. [8t^c)-[8t^e). Note that the cross sections for gg — )■ Z — t- NN are the same as those for 
qq ^ Z ^ E^E^ owing to the same assignment of couplings in Eq. El therefore, they are 
not shown here. 
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FIG. 8: NLO production cross sections for (a) E~^E~ production via an intermediate virtual 
photon and (b) via an intermediate Z boson at the LHC as a function of tue, and (c-e) PDF and 
scale uncertainties of E^E~ production via an intermediate virtual photon. Production via an 
intermediate Z boson exhibits similar uncertainties. 

V. SUMMARY 

We present NLO predictions for the cross sections for new heavy quark and lepton pro- 
duction at three Large Hadron Collider energies. We treat the cases of pair production of 
heavy quarks via strong interactions, single heavy quark production via electroweak interac- 
tions, and the production of exotic heavy leptons. Theoretical uncertainties associated with 
the choice of the renormalization scale and the parton distribution functions are specified. 
We derive a simple and useful parametrization of our results which should assist in experi- 
mental simulations and in phenomenological studies of new physics models containing new 
heavy quarks and/or leptons. 
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under Contract No. DE-AC02-06CH11357. Q. H. C. is supported in part by the Argonne 
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TABLE IV: The fitting parameters {A, B, C) in the parametric formula [Eq. [T] for the NLO total 
cross sections (pb) for heavy lepton production at the LHC. The central values are given for the 
scale choice = fip = fJ-o = fnE- The superscript and subscript denote the deviations of the 
coefficients from the central value when the scale is (superscript) or 2^o (subscript). We use 
the CTEQ6.6M PDFs. 



Process 



Parameter 



^= 7TeV 



lOTeV 



^= 14TeV 



E+N 



A 
B 
C 



0.71010 



+0.00073 
+0.00943 



0.77778 



-0.00278 
+0.00406 



-2.905051q;q5227 



-9.41140: 



-0.05659 
-0.09655 



ifi7i 7+0.02513 
-0.lU(i(_QQ2Q97 

-6.66323i°:S 



0.81919 



-0.00250 
+0.00335 



oil 7/17+0.02228 
-0.il /^/_Q^Q2]_ 79 



-5.20543: 



-0.02221 
-0.02455 



E-N 



A 
B 
C 



0.99258 



-0.00287 
+0.00593 



0.98794 



-0.00431 
+0.00560 



-5.03622 



+0.02662 
-0.03759 



o r;qi c;/|+0.05327 



-4 746S1+0-03171 

4. ^41D01_Q Q33gQ 

-6 4Sq70+°-°2594 

D.^OyrU_Q^Q2g74 



0.94199 



-4.18717" 



-0.00381 
+0.00486 

-0.02669 
-0.02725 



-5 30307+0.01772 
'^•"-'^'-'^'-0.01906 



E+E- 



A 
B 
C 



0.90798 



-0.00402 
+0.00298 



0.89357 



-0.00376 
+0.00524 



^ 7/IS1S+0-03299 
-5.('48i8_oo2533 

S 1 8'^RS+0-04472 



-5 46629+°-°3°^° 

O.^DDZy_g Q3321 

fi 900^:7+0. 02602 
-D.ZOOO^_g Q25gg 



O.86O8O: 

-5.01876" 



-0.00328 
r0.00425 

-0.02561 
-0.02570 



^ 1009^+0-01799 



qq^ Z ^ E+E- 



A 
B 
C 



0.94114 



-0.00508 
+0.00760 



-3.56798 



+0.03728 
-0.04573 



0.90624;°:^ 
-3.22104l°:S 



-8.45564l°;S 



-6.44243 



+0.02677 
-0.02707 



0.88538 



-0.00339 
h0.00434 



O 0071 f;+0.02556 
-Z.5Z;iD_Q^02557 



-5.28945: 



-0.01899 
-0.02084 
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Note_added: While finalizing the write-up of this work, we became aware of a recent 
in which the t channel production of single t' is calculated at the LHC. Our 



paper 



results agree with those in Ref. 
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Appendix A: Standard model parameters 



For our numerical evaluations, we choose the following set of SM input parameters: 
Gf = 1.16637 X lO^^GeV"^ a = 1/137.0359895, 



Following Ref. 



mz = 91.1875 GeV, 
mt = 173.1 GeV, 



Osimz) = 0.1186, 
sin^ = 0.2314. 



we derive the W boson mass as mw = 80.385 GeV. Thus, the square of 



the weak gauge coupling is = 4:\/2m^^GF- 



Appendix B: Numerical fit of the heavy quark pair production cross section 

In this appendix we use heavy quark pair production via the QCD interaction as a sample 
illustration of how one may extend the region of applicability of a fit to a wider mass range 
simply by increasing the number of terms in the polynomial expression. We treat heavy 
quark production with masses within the range 250 GeV and 700 GeV. We choose the 
following formula to fit the cross section: 



log 



a(mQ,/i) 



^+- + C + Dx + Ex\ (Bl) 



X 



pb 

where x = mg/TeV. The fitting functions at the LHC at three energies are given as follows: 

7TeV : 1.33969x^ - 10.4918x + 1.36903 + 1.32168/x - 0.0684983/x^ 
10 TeV : 1.87135a;2 - 9.82576a; + 2.34236 + 1.28886/a; - 0.0682711/x^ 
14 TeV : 2.16753^^ - 9.46830x + 3.33295 + 1.22263/x - 0.0647251/x^ (B2) 

where the renormalization and factorization scales are chosen a.s fiR = fip = fJ'O = ^q- 
Table IVl displays the exact and fitted values of the NLO cross section (pb). The numerical 
values of the exact and fitted cross sections agree very well. 



Appendix C: Influence of parton luminosities on the heavy fermion mass depen- 
dence 

We comment in this appendix on the origins of the remarkably simple analytic expres- 
sion Eq. [1] and on its range of applicability in the heavy fermion mass. The expression 
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TABLE V: The exact and fitted values of the NLO cross section (pb) for heavy quark pair pro- 
duction at the LHC at three cm. energies, where the scale is chosen as = fip = fiQ = mq. 



(GeV) 


7TeV 
Exact Fit 


lOTeV 
Exact Fit 


14TeV 
Exact Fit 


250 


20.51 20.50 


58.34 58.32 


142.10 142.06 


300 


7.289 7.291 


22.21 22.21 


57.03 57.04 


350 


2.940 2.939 


14.38 14.38 


25.78 25.78 


400 


1.301 1.301 


9.565 9.561 


12.74 12.74 


450 


0.617 0.618 


2.285 2.286 


6.742 6.743 


500 


0.310 0.310 


1.224 1.224 


3.771 3.770 


550 


0.162 0.162 


0.685 0.685 


2.203 2.204 


600 


0.088 0.088 


0.398 0.398 


1.336 1.337 


650 


0.049 0.049 


0.239 0.239 


0.837 0.837 


700 


0.028 0.028 


0.147 0.147 


0.539 0.539 



works well for strong pair production of heavy quarks, as well as for electroweak single 
heavy fermion production and electroweak pair production of heavy fermions. This appar- 
ent process independence suggests that the heavy fermion mass dependence of the parton 



luminosities 



40- 



- 42| is playing a dominant role. 
The quark-antiquark and gluon-gluon contributions to the cross section for heavy quark 
pair production at hadron colliders can be written as 



E 



dr 



Xn 



fq/P (a^a,/io) fq/I 



Xn 



S(^qq^QQ\S 



)] 



/ '^fa/P (a^a, /io) fg/p /io) T [saqq^QQ{s)] 



(CI) 



where a is the operative parton-level cross section, r^j^ = 4mQ/s, and s = ts with s 
being the square of the cm. energy of LHC. In order to represent a process independent 
situation, we approximate the (dimensionless) square brackets by one. As a result, the only 
remaining effect of the heavy quarks appears in the lower limit of integration Tmin- We 
choose fiR = fip = fng and use the CTEQ6.6M average PDF set. 

Using Eq. \C1\ we generate "cross sections" as a function of rriQ, and we then fit the 
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FIG. 9: Fit to the trial scattering cross sections in the range of 0.5 < uiq < 2.0 TeV: (a) for the 
quark-antiquark annihilation channel, and (b) for the gluon fusion channel. 

results with Eq. [1] The agreement is excellent for the qq and gg cases over the range of 
heavy quark masses shown in Fig. [9l 0.5 < mg < 2.0 TeV. The square symbols in the figure 
are the predicted cross sections, and the solid curves represents the fits. 

If we try to extend our study to smaller masses than shown in Fig. [9] we find that the 
simple expression Eq. [1] is no longer adequate, as shown in Fig. [TOTa). When Tmm is too 
small, we sample parton densities in the region of small Bjorken x where they grow rapidly 
with decreasing x. More terms in the polynomial are needed for a better fit. We can fit the 
gluon fusion cross section well with the addition of two more terms; see Fig. [TOT b) 
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FIG. 10: Fit to the trial scattering cross section for the gg channel extended into the region of 
small mq. 
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